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Abstract We discuss the new integrable boundary conditions for the 0{N) non- 
hnear a model and related solutions of the boundary Yang-Baxter equa- 
tion, which were presented in our previous paper hep-th/0108039. 
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1. Introduction 

Two-dimensional nonlinear a (nlu) models have been the subject of 
intense study during the past few years, since they may be used as 
toy models for the study of higher dimensional non-abelian gauge theo- 
ries (Yang-Mills), they arise in several condensed matter and statistical 
mechanics problems, and there are powerful mathematical methods in 
2d that allow one to have a deeper understanding of their structure. 
On top of that, they display a host of theoretical phenomena, such as 
asymptotic- freedom, dynamical mass generation, and 1/iV-expansions. 

A natural generalization, both from the theoretical and experimental 
point of view, of a given 2d integrable model, is to consider it on the half- 
line [1]. This type of reduction arises in several problems, for example, 
when considering the radial part of the Schroedinger equation for a radial 
potential, in the study of quantum impurities, such as the Kondo effect, 
and in open string theory. 

In this note we summarize the results obtained in [2], where we have 
found new integrable boundary conditions and related solutions of the 
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boundary Yang-Baxter equation (bYBe) for the 0{N) n\a model. The 
following discussion is informal and intended to a general audience, we 
refer to [2] for a more complete discussion. 

2. The 0{N) Nonlinear Sigma Model 

The lagrangian of the 0{N) nla model [3] is 

Cnla = ^{d^f (1) 

where the field n = (ni, n2, . . . , njv) satisfies = 1, and is a coupling 
constant. This constraint may be introduced in the lagrangian through 
a lagrangian multiplier A(a;), 

Cnia = ^{dnf + m{i^^-^) (2) 

The equation of motion for the field n is then easily derived from (2) 
and, in light-cone coordinates x± = (xq ± xi)/2, it reads 

5+a_n = -n • a_n) (3) 

The exact S'-matrix for te 0{N) nla model was found by Zamolodchikov 
and Zamolodchikov in [4] . For a review on exact S'-matrices see [5] . 

Since the 0{N) nlo" model is classically conformally invariant (no di- 
mensionful quantities appear in the classical lagrangian (2)) the trace 
of the energy-momentum tensor, (= T_+), vanishes. Upon quan- 
tization, conformal invariance is broken due to the introduction of an 
energy scale, such as an energy cut-off', for example. In any case we will 
see later that the classical conformal invariance is crucial in establish- 
ing the quantum integrability of this model. In light-cone coordinates 
the only non-vanisihing components of the energy momentum tensor are 
r+_|_ = 9+n • 9+n and T = d-H ■ d-H. This means that energy- 
momentum conservation reads 

d-T++ = (4) 

and an analogue equation with + — . Classically it is easy to sec that 
(4) implies 5_(r+_|_)" = for any integer n. Upon quantization this is 
no longer true since we have now a composite operator and one must 
be careful in defining the product of operators at the same space-time 
point. This means that the right-hand side of (4) will no longer be zero, 
or in other words, we have a quantum anomaly. Even though it is a hard 
problem to compute the exact form, with all numerical coefficients, of the 
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anomaly in (4) , we can greatly fix its form, by using dimensional analysis 
and group theory. Let us look at what happens in the case n = 2. In 
this case the rhs of (4) has dimension 5, Lorentz weight 3 and it is a 
scalar under the global 0{N). This means that, whatever the anomaly 
is, it is a local operator that satisfies these three requirements. One can 
then proceed and make a list of the possible operators that contribute to 
the anomaly. The rhs of (4) will be, finally, a linear combination of these 
operators, the hard work being to have to compute the coefficients of 
this linear combination. In the case of 5_(T+_|_)^ = a wonderful thing 
happens: all the operators that can contribute to the anomaly can be 
rewritten as a total derivative, with the help of the equation of motion 
and the constraint n • n = 1. This is the so-called Goldschmidt-Witten 
argument [6] (see also [3], for an earlier version). For the complete list 
of these operators, we refer to [6]. By using this quantum conservation 
law the integrability of the 0{N) nla model is established. 

When considering a field theory on the half-line, bulk conservation 
laws may be broken, as it is clear in the case of linear momentum con- 
servation (the boundary breaks translation symmetry). Therefore, one 
needs to impose suitable boundary conditions that will ensure that a 
given bulk conservation law will still hold after the introduction of the 
boundary. In equations, if we have a bulk conservation law of spin 
(Lorentz weight) s 

5+ Ji^+^) = and d-J^^+'^ = d+R^r'^ (5) 

it follows that the charges 

Q±= r^dxi{jt'^-R^4-'^) (6) 

J —oo 

are conserved. After the introduction of a boundary, only (possibly) a 
combination of these charges will be conserved. The charge is 

Q = r dx, {4+'^ - + Ji^+') - i?J-')) + m (7) 

J —oo 

where S(i) is a local operator satisfying 

^(.+ 1) _ ji^S+l) ^ j^(s-l) _ 

This is precisely where we have to impose the boundary conditions in 
order to have a non-trivial conserved charge in the presence of a bound- 
ary. 
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In [7] we have shown that if we impose Neumann boundary conditions, 
dini\x=o = 0, to field components, and Dirichlet boundary conditions, 
(?onj|x— = 0, to the remaining N — k components, the condition (8) is 
satisfied and we have integrable boundary conditions. Note that there is 
no free parameter (coupfing constant) in this case. These are diagonal 
boundary conditions, in the sense that scattering off the boundary does 
not change the 0{N) index of the incoming particle. These boundary 
conditions break the bulk symmetry at the boundary to 0{k) x 0{N—k). 
It can be shown, by using the bYBe, that these are the only possible 
diagonal integrable boundary conditions. Therefore, if we are looking 
for new integrable boundary conditions with free parameters, we have to 
necessarily look for non-diagonal boundary conditions, that is, boundary 
conditions that change the 0{N) index (flavor) of the incoming particle 
through scattering off the boundary. 

In the following we will take a slightly different point of view from [2] , 
where we used a simple two free bosons model coupled at the boundary 
as a guide to the possible boundary conditions for the 0{N) nla model. 

In [2] we have found new integrable boundary conditions that break 
the bulk 0{N) symmetry to 0(2) x 0{N - 2) at the boundary, and 
which depend on one free-paranictcr g. The reason for this symmetry 
at the boundary is the following. Free boundary conditions (Neumann) 
have 0{N) boundary symmetry. The diagonal boundary conditions we 
found previously, break the boundary symmetry to 0{k) x 0{N — k). 
We are looking now for non-diagonal boundary conditions with a free 
parameter, and we may assume that in certain limiting cases, such as 
taking the coupling constant to or oo, we should reduce to a diagonal 
case. Therefore, if we insist that once the boundary symmetry is broken 
wc do not have any point in the integrable flow where it is enhanced, 
we should look for non-diagonal boundary conditions that are 0{k) x 
0{N — k) symmetric. By writing the generic boundary condition as 

dini\x=o = Mij donj\^=o (9) 

where the indices i and j run through a subset of {1,2,... , A^} (the 
''nondiagonal subset"' , which could be taken to be the first k indices, for 
example). Under an orthogonal transformation the fields transform as 
ni ^ ha = OaiTii, where O is a kxk orthogonal matrix. This means that 
Mij —I- Mab = OaiMijOhj, and if we require the boundary conditions to 
be 0{k) x 0{N — k) symmetric, we should have OMO^ = M. The only 
case where we can impose this condition for a non-diagonal matrix if 
when k = 2, since 0(2) is abelian. This fixes the matrix M to be of the 
form M = gil + i §2(^2, where / is the identity matrix and 02 is a Pauli 
matrix. By inspecting the spin-4 Goldschimdt-Witten charge described 
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above we see that if we take 51 = and §2 = 9 arbitrary, the following 
boundary condition is integrable: 

dini\x=o = g don2\x=o and din2\x=o = -g doni\a;=o (10) 

where we picked the first two components of the n field without any loss 
of generality. The remaining field components satisfy Dirichlet bound- 
ary condition ^. In a different form, this boundary condition had been 
studied by Corrigan and Sheng at the classical level in [8], for the 0(3) 
nlcr model. 

The non-diagonal boundary conditions in (10) can be derived from 
the boundary lagrangian Cj, = ^MijUifij, which shows that Mij should 
be anti-symmetric. 

By taking 5 — > we have diagonal boundary conditions, where 2 field 
components satisfy Neumann and the remaining Dirichlet, and by taking 
g ^ 00 we recover a diagonal case again, with all field components sat- 
isfying Dirichlet boundary conditions. Therefore wc have an integrable 
flow between diagonal boundary conditions, from 0{N), corresponding 
to = 00, to 0(2) X 0{N — 2), corresponding to g = 0. 

3. The Reflection Matrix 

When one tries to find an exact S'-matrix for a given integrable field 
theory, the use of the bulk symmetries plays a crucial role, making it 
much easier to solve the Yang-Baxter equation. This is why we had to 
understand the symmetry of the boundary conditions before we could 
go on and try to solve the bYBc. 

For the purely diagonal case, the solutions of the bYBe have been 
found in [9]. They are block diagonal, 0{k) x 0{N — k) symmetric, 
with diagonal elements {Ri{9), . . . , R2{0) . . .), the first k elements corre- 
sponding to Neumann, and the remaining N — k to Dirichlet. The bYBe 
fixes the ratio Ri{0)/R2{6) to be 

Ri{e) c-e 



R2ie) c + i 



(11) 



with c = — if ^n^2^ - Note that there is an interesting duality by taking 
k ^ N — k, which takes c —>■ — c, and therefore Ri{9) R2{Q)- 

For the boundary conditions (10), we start with the following ansatz 

/ A{e) B{0) ••• \ 

-Bie) A{e) 

R= Roie) ■■■ . (12) 

Roie) ■■■ 
■ ■■ J 
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This means that the first two particles can scatter onto each other with 
ampUtude ±B{6), or onto themselves with amplitude A{9). The di- 
agonal elements correspond to the particles scattering diagonally with 
Dirichlet boundary conditions, with amplitude Ro{9). Thinking in terms 
of the boundary lagrangian for the non-diagonally scattering particles, 
we see that the off-diagonal amplitudes should have opposite signs. 

We can use the bYBe now, in order to compute the functions A{9), 
B{6), and Ro{6). In the following we will quote the differential equations 
for X{9) = A{e)/Ro{e) and Y{e) = Ble)/Ro{e), obtained from the 
bYBe, by taking the limit where the two rapidities are equal. 

The process \Ai{ei)Ai{e2)) \Ai{-0i)Ai{-92)) ^, where i is any of 
the diagonally scattering particles, gives 

The process \AMMO2)) ^ \Ai{-eM2{-e2)) gives 

These two equations can be easily solved by the introduction of the 
auxiliary functions Z±{9) = X{6) it iY{6). We obtain 

1/c-e c'-0\ , 1 /c-e cZ-i 

(15) 

where c and c' are constants to be determined. Since we have only one 
free parameter at the boundary, we should find one equation relating c 
and c'. This is acomplished by the bYBe corresponding to the process 
|^i(6'i)-4i(6'2)) |^i(-6'i)^2(-6'2)), from which we obtain 

N - 4 

c + c = -^TT^^— ^ (16) 

We have verified that with this constraint, all the other bYBe's are 
satisfied. Once the ratios X{9) and Y{9) have been fixed, all that is left 
to do is to compute the overall factor for the reflection matrix, which 
can be done with the use of boundary unitarity and boundary crossing- 
symmetry, and a minimality hypothesis for the pole structure of the 
reflection matrix. We refer to [2] for the explicit results. 

Note that if c = c' the off-diagonal amplitudes ±Y(6) vanish, and 
we recover a diagonal scattering problem. The other instance where 
Y(9) vanishes is when |c|, |c'| —i- 00. In the first case the ratio X{0) = 
with c = — ^E|) which corresponds precisely to the case where 
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the first two field components satisfy Neumann boundary conditions 
and the remaining N — 2 Dirichlet. This is the reason why we chose 
the diagonally scattering field components to satisfy Dirichlet boundary 
conditions. In this case the solution for the reflection matrix is 0{2) x 
0{N — 2) symmetric. By looking at the explicit form of the boundary 
conditions (10), wc see that this corresponds to ^ = 0. In the second 
case X{6) = 1, which means that the reflection matrix is proportional to 
the identity, and therefore 0{N) symmetric. This corresponds to taking 
g ^ oo, and therefore, to all components satisfying Dirichlet boundary 
conditions. 

We can introduce the following conveninent parametrization: c = 
— + ^(q) and c' = — 'jf ~ ^(5)) where ^{g) is an unknown 
function of the boundary coupling constant. The two cases described in 
the preceding paragraph correspond to ^(0) = and $^{g 00) 00. 
This establishes an integrable flow between different diagonal boundary 
conditions. 

One could be tempted at trying a generalization of the ansatz (12), 
with more than one non-diagonal block, corresponding to more than one 
pair of particles being coupled at the boundary. By using the bYBe it 
can be shown that there are no solutions of this type [2] . 

In [2] we found other solutions to the bYBe for the 0{2N) nla model, 
but were not able to link them to any boundary conditions. Another 
special case is the 0(2) nlcr model. Naively one could be lead to think 
that the 0(2) nla model is equivalent to a masslcss free boson, through 
a mapping (ni,n2) (cos(^), sin(^)), but this is not the case, and after 
a more careful analysis, it can be shown that the 0(2) nlcr model is 
equivalent to the sine-Gordon model at = Svr, which describes the 
Kosterlitz-Thouless point of the classical XY model. The solution we 
found depends on three parameters, instead of two as in the boundary 
conditions found by Ghoshal and Zamolodchikov in [1]. The resolution 
of this discrepancy is that we are looking at the sine-Gordon model at 
a special value of the coupling constant, and as already remarked in [1], 
at these special points there are more solutions than the ones found for 
the general case. 

4. conclusions 

We have found new integrable boundary conditions for the 0{N) nlcr 
model, which depend on one free parameter g. These bounday conditions 
break the bulk 0{N) symmetry to 0(2) x 0{N — 2), and by taking 
the limits g ^ and 5 ^ 00 we recover diagonal solutions studied 
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previously. This establishes an integrable flow between two different 
sets of boundary conditions. 

Recently Mackay and Short [10] have studied the principal chiral 
model with a bounday, and found an interesting relationship between 
their boundary conditions and the theory of symmetric spaces. Their 
solutions, though, are quite different from ours, and some work should 
be done in trying to clarify their relationship. 

As natural follow-up problems, one should try to find explicitly the 
function ^{g) in the reflection matrices, and to study the boundary ther- 
modynamic Bethe ansatz equations. 

An interesting direction to pursue would be to extend these results 
to the SO{N) Gross-Neveu (GN) model. Since the S'-matrix for the 
elementary fcrmions of the GN model is equivalent to the one for the 
0{N) nla model, up to a CDD factor, we certainly can find solutions of 
the bYBe of the form (12) for the GN model too. 
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Notes 

1. If a quantity tp transforms a,s (j> ^ exp{pa)<p under a Lorentz boost, that is, the rapidity 
variable 6 — > 6 + «, we say that (p has Lorentz weight p. 

2. This choice of boundary condition for the remaining field components will become clear 
when we discuss the boundary Yang-Baxter equation. 

3. The {Ai{6)} are the usual Faddeev-Zamolodchikov operators. 
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